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A time-domain approach is proposed for the propagation of ultrashort electromagnetic wave pack- 
ets beyond the paraxial and the slowly-varying-envelope approximations. An analytical method 
based on perturbation theory is used to solve the wave equation in free space without resorting 
to Fourier transforms. An exact solution is obtained in terms of successive temporal and spatial 
derivatives of a monochromatic paraxial beam. The special case of a radially polarized transverse 
magnetic wave packet is discussed. 



State-of-the-art mode-locked laser oscillators can gen- 
erate ultrafast pulses whose duration can be as short as 
a few optical cycles 0. Such pulses - with a finite spa- 
tiotemporal extent and thus carrying a finite quantity of 
energy - can in principle be isolated, amplified, tempo- 
rally compressed, and spatially focused to reach extreme 
field intensities within a volume that is of the order of the 
third power of the laser wavelength (A 3 ) 1]. In this sit- 
uation the commonly used paraxial approximation and 
slowly-varying-envelope approximation (SVEA) do not 
apply directly and must be corrected. 

In 1975, Lax et al. proposed a method to express exact 
solutions of Maxwell's equations in terms of an infinite sc- 
ries of corrections applied to a paraxial electric field 0- 
So far, this approach has led numerous studies dealing 
with non-paraxial optics. Unfortunately, as it is also the 
case for the electromagnetic beam theory developed by 
Davies [j| , the so-called "method of Lax" only applies to 
monochromatic beams, or CW laser signals. With the 
outstanding developments in femtosecond and attosec- 
ond technologies, more appropriate tools are required to 
describe the propagation of non-paraxial ultrashort light 
pulses, sometimes called light bullets. Recently, Lu et 
al. Q have proposed a method to describe the free-space 
propagation of vectorial non-paraxial ultrashort pulsed 
beams in the frequency domain. In this letter, we present 
an alternative approach to this problem entirely in the 
time domain, i.e. without using Fourier transforms. 

The spatiotemporal propagation of an electromagnetic 
wave packet is governed by Maxwell's equations. Solving 
these equations in free space proceeds in three steps : 
1. the transverse electric field vector Ex is found from 
the wave equation, 2. the longitudinal electric field E z 
is obtained from the divergence equation (V • E = 0), 
and 3. the magnetic flux vector B is calculated from the 
differential form of Faraday's law (V x E = —9(B) or 
the differential form of Ampere's law (VxB = c -2 <9 t E), 
where c is the speed of light in vacuo, and dt = d/dt. 
In the case of a paraxial wave packet whose duration 
is much longer than the optical period, this sequence of 
operations is performed with ease. We now describe each 
of these three steps in detail beyond the range of validity 
of the paraxial approximation and of the SVEA. 

In free space, the vectorial wave equation for the trans- 



verse electric field vector Ex - perpendicular to the time- 
averaged Poynting vector - is given by the following wave 
equation 
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The electric field vector is defined as follows in phasor 
notation (see also Refs. G] and0) : 
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where Ex is the complex envelope of the transverse elec- 
tric field vector, E z is the complex envelope of the longi- 
tudinal electric field, a z is a longitudinal unit vector ori- 
ented along the pulse propagation, j = V~-l ; ujq = hoc, 
fc = 27rAo, and Ao are the central angular frequency, the 
central wave number, and the central wavelength of the 
wave packet spectrum, respectively. Combining Eqs. 
and (J2J, one obtains the following equation for the com- 
plex envelope of the field : 
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where Vx is a differential operator acting in the trans- 
verse plane. 

The use of the retarded time t' = t — z/c instead of 
the (normal) time t is sometimes preferred when dealing 
with the propagation of electromagnetic signals of finite 
spatiotemporal extent p|. More specifically, it separates 
the transverse envelope (diffraction) from the temporal 
envelope (pulse shape). After some manipulations, the 
result given at Eq. (0 can be expressed in terms of the 
retarded variables t' = t — z/c and z' = z as follows : 

VlEx - 2jk d z ,E ± + 9 2 ,Ex - - d z ,d t ,E x = 0. (4) 
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Alternatively, Eq. can be written in a more compact 
form that reads : 



ViEx-2jfco9^[l-e]Ex = 0, 
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with O — j (oj 1 dt> — (2fco) 1 9 2 '). For paraxial and 
quasi-monochromatic optical beams, the term 0Ex is 
vanishingly small compared to the complex field Ej 



it- 



self, 
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Ex- In that case, Eq. (J3J leads 

to the paraxial wave equation 0. On the other hand, 
when the transverse (or longitudinal) dimension of the 
wave packet is of the order of the wavelength, the term 
©Ex can no longer be neglected. However, the modulus 
of this contribution, i.e. |0Ej_|, remains small even for 
strongly focused single-cycle light pulses 0, IE IE • As 
a consequence, we can consider Eq. (j5J) to be a perturbed 
paraxial wave equation. 

We observe that the partial differential operator © 
cannot be reduced to a simpler operator, proportional 
only to spatial or temporal derivatives, i.e. O cx d z or 
cx <9™. This fact clearly indicates that spatial and 
temporal effects cannot be decoupled. Consequently, a 
solution of Eq. is found by expanding the transverse 
electric field vector Ex as a power series of 0. In a com- 
pact notation it reads : 
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z' to yield : 
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A general equation giving the magnetic flux vector B 
of an arbitrary ultrashort wave packet can be deduced 
from Maxwell's equations but cannot always be reduced 
to a simple expression, as it is the case for monochro- 
matic transverse electromagnetic (TEM) waves (see Eq. 
(7.11) of Ref.i|. However, for a given distribution and 
polarization of the transverse electric field vector corre- 
sponds a spatiotemporal arrangement of the fields that is 
unique and guarantees the stability of the wave packet. 
As an example, let us consider the case of a TMoz wave 
packet whose transverse electric field is radially polar- 
ized, i.e. Ex = E r a r and Eg = 0, where r and 9 are polar 
coordinates in the plane perpendicular the propagation 
axis (the z-axis), and I = 1, 2, 3, . . . The intensity profile 
of this particular family of beams is characterized by I 
concentric and angularly symmetric (dgE± = 0) bright 
rings 0. From Ampere's law (VxB = c~ 2 <9tE), we ob- 
tain the two following equations (B r = and B z = 0) : 



By equating terms with the same power of in Eq. JSJ, 
the two following equations are then found : 

Vi* ( L 0) -2. ? fc 5 2 '* ( L ° ) = 0, (7a) 

vi#i n) - 2jk d z ,v<? ) + 2jk Q d z ^- 1] = 0, (7b) 

solving the perturbed paraxial wave equation in such a 
way that the nth contribution is obtained recursively 
from the order (n — 1), acting like a source term in the 



paraxial wave equation (\I/ 
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A solution of the recursion relation given at Eq. J7 
is written as follows (see Porrasjg) : 
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thus giving in terms of successive derivatives of 

the zero-order field , a solution of the paraxial wave 
equation. Recalling Eq. ©, the envelope of the trans- 
verse electric field vector reads : 
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For any given transverse electric field vector Ex, the 
longitudinal electric field E z can be obtained from the di- 
vergence equation. With the use of the complex notation 
defined at Eq. © , V • E = can be formally inverted 01 
and expressed in terms of the retarded variables t' and 
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If we combine Eqs. <|lla|) and (|llb|) together, we obtain 
the divergence equation that we solved earlier at Eq. llUfl . 
Thus, the solution of Eq. (|llafl assures that Eq. (|llbl) 
is also respected. If we express the fields in complex 
notation (and then in terms of the retarded variables), 
the formal inversion of Eq. l|llaf) yields : 
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In the past, many investigations have been devoted to 
the study of the paraxial wave equation and its solutions 
(see for example Ref. ITTI) . More specifically, we make 
reference to the wide family of the Gauss-Hermite and 
Gauss-Laguerre beam modes, considered to be a good 
representation of the transverse electromagnetic modes 
of optical resonators 0, 0|- According to the formal- 
ism we propose, these solutions, or any solution obtained 
under the paraxial beam optics theory, can be extended 
to the case of tightly focused (non-paraxial) ultrashort 
(few-cycle) wave packets using the general formulation 
we developed. We emphasize the fact that we obtained 
all the spatiotemporal corrections in the time-domain, al- 
ternatively to the Fourier transform approach 4J. How- 
ever, the method described here has the feature that it 
displays directly the corrections of the SVEA on the four- 
dimensional structure of the wave packet, as needed to 
correctly predict ultrafast interactions with matter. For 
instance, dispersion and nonlinear optical effects could 
be embedded in this formalism [(J. 
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This letter should be considered as an extension - a 
vectorial and non-paraxial generalization - of a previ- 
ous work published by Porras In fact, it is possi- 
ble to show that Eq. © simplifies to the equation for 
paraxial ultrashort pulsed beams when the spatial varia- 
tions of the field due to diffraction are negligible com- 
pared to the variations of the temporal envelope, i.e. 
when 6 ~ juj^ 1 dt' (see Eq. (8) of Ref. |3). Besides, it 
has been demonstrated that the perturbation method we 
used leads to the same results that are obtained from 
integral methods Q- We thus observe that the fields 
given at Eqs. 10 and (|1U[) are equivalent to the Eqs. (25) 
and (28) of the paper of Lu et al. Q. 

The irreducibility of the expansion parameter in- 
dicates that in free space, the spatial and temporal en- 
velopes of a non-paraxial ultrashort wave packet are non- 
separable. This coupling ensures the wave packet's sta- 
bility under tight focusing and reveals itself as a reorga- 
nization of the structure of the wave packet, noticeable as 
modifications of the carrier wave and the spatiotemporal 
envelope @. 

For a gaussian beam (with a gaussian pulse 
shape), the spatial and temporal contributions, i.e. 
(2/c )~ n <9™, and w^""9™Ex, are respectively propor- 
tional to (2k z R )~ n = (fc r )~ 2n and (cj Ai)~ n , where 
Z R = fco^o/^ is the Rayleigh distance, r$ is the wave 
packet transverse extent at focus (beam waist), and At 
is the pulse duration. The expansion of 9™ yields also 



crossed-terms with various powers of k^r® and uj^At. 
However, for a single-cycle pulse (At = 2tt/u>q) focused 
down to one wavelength (tq = 27r/fco)j the most im- 
portant contribution comes from the temporal envelope 
and decreases as (2n)~ n . Consequently, it is observed 
that a truncated series (n = 0, 1, 2) reproduces the exact 
field with a good accuracy (see Porras 8J). Nonetheless, 
for slightly focused and long pulses (ro >> 2ir/ko and 
At >> 27r/a;o), the formalism reduces to the expressions 
for paraxial and monochromatic beams 0, lllj . 

The full vectorial treatment we have presented here 
also leads to an exact evaluation of the associated longi- 
tudinal electric and magnetic fields. These field compo- 
nents are not usually dealt with in optics ; however, they 
must be taken into consideration for the study of rela- 
tivistic effects in laser-matter interactions and for the in- 
vestigation of electron acceleration in intense laser fields. 
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